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THE SPEED AND ORIENTATION OF THE PARSEC-SCALE JET IN 3C 279 

B. Glenn Piner 1 , Stephen C. Unwin 2 , Ann E. Wehrle 3 , Alma C. Zook 4 , C. Megan Urry 5 , 

& Diane M. Gilmore 6 

ABSTRACT 

A high degree of relativistic beaming is inferred for the jets of blazars based on several lines of evidence, 
but the intrinsic speed and angle of the jet to the line-of-sight for individual sources are difficult to 
measure. We have calculated inverse- Compton Doppler factors for 3C 279 using the collection of VLBI 
data (including high-resolution space VLBI data at low frequencies) recently published by us (as Wehrle 
et al. and Piner et al.), and the collection of multiwavelength spectra recently published by Hartman 
et al. From the Doppler factor and superluminal apparent speed, we then calculate the Lorentz factor 
and angle to the line-of-sight of the parsec-scale relativistic jet. We follow the method previously used 
by Unwin et al. for 3C 345 to model the jet components as homogeneous spheres and the VLBI core as 
an unresolved inhomogeneous conical jet, using KonigPs formalism. 

The conical-jet model can be made to match both the observed X-ray emission and the VLBI properties 
of the core with a suitable choice of Doppler factor, implying the core makes a significant contribution 
to the X-ray emission, in contrast to the situation for 3C 345, where the jet components dominated the 
X-ray emission. The parameters of the Konigl models indicate the jet is particle dominated at the radii 
that produce significant emission (from ~ 5 to 20 pc from the apex of the jet for most models), and is 
not in equipartition. At the inner radius of the Konigl jet the magnetic field is of order 0.1 G and the 
relativistic-particle number density is of order 10 cm -3 . The kinetic energy flux in the jet is of order 
10 46 (1 + k) ergs sec -1 , where k is the ratio of proton to electron energy, which implies a mass accretion 
rate of order 0.1(1 + k)/rjMQ yr -1 , where rj is the efficiency of conversion of mass to kinetic energy. 

When all components are included in the calculation, then on average the core produces about half 
of the X-rays, with the other half being split between the long-lived component C4 and the brightest 
inner-jet component. We calculate an average speed and angle to the line-of-sight for the region of the 
jet interior to 1 mas of v = 0.992c (7 = 8) and 9 = 4°, and an average speed and angle to the line-of-sight 
for C4 (at r « 3 mas) of v = 0.997c (7 = 13) and 8 — 2°. These values imply average Doppler factors of 
S = 12 for the inner jet, and 6 — 21 for C4. 

Subject headings: quasars: individual: (3C279) — galaxies: jets — galaxies: active — radiation 
mechanisms: non-thermal — radio continuum: galaxies 



1. INTRODUCTION 

The quasar 3C 279 (z=0.536) has been one of the most 
intensively studied quasars for several reasons. It was the 
first radio source observed to exhibit the phenomenon of 
apparent superluminal motion (Knight et al. 1971; Whit- 
ney et al. 1971; Cohen et al. 1971), prompting continued 
study with VLBI through the 1970's and 1980's (Cotton 
et al. 1979; Unwin et al. 1989; Carrara et al. 1993) 
and 1990's (Wehrle et al. 2001; hereafter Paper I). A 
bright 7-ray flare was observed from 3C 279 in 1991 by 
the EGRET instrument shortly after the launch of the 
Compton Gamma Ray Observatory (Hartman et al. 1992), 
leading to a great deal of multiwavelength coverage during 
the following nine years. 3C 279 was one of the brightest 
7-ray blazars detected during the lifetime of the EGRET 
instrument (Hartman et al. 1999). It is one of the proto- 
types for the class of luminous "red" blazars (Sambruna 
2000). 



These studies have produced a large amount of data on 
3C 279. A number of nearly simultaneous multiwavelength 
spectra are available that show 3C 279 at various levels 
of activity. Data previously published by Maraschi et al. 
(1994), Hartman et al. (1996), and Wehrle et al. (1998) 
are compiled by Hartman et al. (2001a) (hereafter H01). 
Strong variability on timescales of a day or less has been 
observed in optical through 7-ray bands (e.g., Wehrle et 
al. 1998; Hartman et al. 2001b). In the radio regime, the 
variability timescale is longer, and flux density monitoring 
at 4.8, 8.4 and 14.5 GHz, complete with polarization data, 
has been obtained at the University of Michigan Radio Ob- 
servatory (Aller et al. 1985). Monitoring at 22 and 37 GHz 
has been done at Metsahovi Observatory (Terasranta et al. 
1992, 1998). We have recently published a compendium 
of six years of VLBI images of 3C 279 at 22 and 43 GHz 
from 1991 to 1997 (Paper I), showing the kinematics of the 
parsec-scale jet. In this paper we combine the VLBI data 
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from Paper I with multiwavelength spectral information 
to calculate the Doppler factor, orientation, and speed of 
3C 279's parsec-scale jet. 

The two-humped overall spectral energy distribution of 
blazars is most naturally explained with a jet of relativis- 
tic electrons emitting a combination of synchrotron radia- 
tion for the radio through optical-uv region, and inverse- 
Compton emission at higher energies (e.g., Wehrle 1999). 
This interpretation of the spectrum leads to derivation 
of physical conditions in the jet that include bulk rela- 
tivistic motion. Evidence for bulk relativistic motion in 
blazars comes from many sources, including the observed 
apparent supcrluminal motions in blazar jets (Vcrmculcn 
& Cohen 1994), the transparency of blazar cores to high- 
energy 7-rays (Dondi & Ghiscllini 1995), rapid flux vari- 
ability (Lahteenmaki & Valtaoja 1999), high VLBI core 
brightness temperatures (Tingay et al. 2001), an excess of 
predicted over observed inverse-Compton emission (Ghis- 
ellini et al. 1993), and arguments invoking equipartition 
or minimum energy requirements (Readhead 1994). 

This bulk relativistic motion can be quantified by the 
ratio of observed to emitted frequency, or Doppler factor, 



7(l-/?cos0)' V ' 

where 9 is the angle to the line-of-sight, (3 = v/c, and 
7 = (1 — Z? 2 ) -1 / 2 is the bulk Lorentz factor. Although ev- 
idence for high values of 5 is strong, calculating S for any 
given source is difficult, as we explain below. Knowledge 
of S is desirable because, together with the apparent super- 
luminal speed, it constrains both the bulk Lorentz factor 
(important for studying jet energetics) and the angle of the 
jet to the line-of-sight (important for unification studies). 
One method of calculating the Doppler factor is to use 
source properties measured from VLBI images and multi- 
wavelength spectra to predict the X-ray flux density that 
should be emitted by the Synchrotron Self-Compton (SSC) 
process. This calculation typically over-predicts the X-ray 
flux density from the source. By assuming the source is 
relativistically beamed with a certain S (see equation (JSJ), 
the conflict between the predicted and the observed X-ray 
flux density can be eliminated (e.g., Marscher 1987; Ghis- 
ellini et al. 1993). If all of the X-rays from the source 
are not due to the SSC process from the component being 
considered, then a lower limit to S is obtained rather than 
a firm value. 

In practice, this method has many problems. It depends 
on the assumed geometry of the emitting region, and it 
depends sensitively on parameters that must be measured 
from multi-frequency VLBI images. Many authors assume 
a homogeneous sphere geometry for the VLBI core (e.g., 
Ghisellini et al. 1993; Mantovani et al. 2000). In re- 
ality, a homogeneous sphere is a poor approximation to 
the VLBI core, because it predicts a sharply-peaked syn- 
chrotron spectrum, and not the flat spectrum over several 
decades in frequency shown by many sources. A more re- 
alistic model for the VLBI core is an inhomogeneous con- 
ical jet, such as the model by Blandford & Konigl (1979) 
and Konigl (1981), which can reproduce the observed mul- 
tiwavelength spectral indices. Estimates for parameters 
such as the synchrotron turnover frequency of VLBI com- 
ponents also bring uncertainty into this calculation, and 



Lahteenmaki, Valtaoja, & Wiik (1999) show that plotting 
SSC Doppler factors derived by various authors from dif- 
ferent data sets for the same sources results in an almost 
pure scatter diagram, most likely due to varying assump- 
tions about input parameters. 

SSC Doppler factors have been calculated in detail per- 
haps only for 3C 345 by Unwin et al. (1994, 1997). They 
modeled 3C 345 as a superposition of an inhomogeneous 
conical jet (using Konigl's 1981 formalism) for the VLBI 
core, and a series of homogeneous spheres for the VLBI 
components or "blobs" . They then used multi-frequency 
VLBI data, quasi-simultaneous X-ray data, and the mul- 
tiwavelength spectrum to constrain the Doppler factor of 
3C 345. In this paper, we apply the procedure used by 
Unwin et al. (1994, 1997) for obtaining inverse-Compton 
Doppler factors to 3C 279. 

Note that our goal in this paper is not to model the 
entire multiwavelength spectrum and short-term variabil- 
ity of 3C 279 using the Konigl model. We use the Konigl 
model as a geometry somewhat more sophisticated than 
a homogeneous sphere to calculate what we hope will be 
an accurate measurement of the SSC Doppler factor, in 
order to constrain the jet orientation and speed. In par- 
ticular, we make no effort to model the 7-ray portion of 
the spectrum, because models including other radiation 
mechanisms, such as that of H01 and Ballo et al. (2002), 
show that this is quite likely due to external Compton 
scattering and not SSC. The model of H01 does indicate 
that the lower energy portion of the inverse-Compton spec- 
trum (the X-rays) is likely to be dominated by SSC emis- 
sion. We also consider only the comparatively "quiescent" 
emission from 3C 279, because the rapid variability during 
flares may originate in components much smaller than the 
Konigl jet, presumably small blobs (like those in the H01 
model) that are moving out through the inhomogeneous 
jet and are superposed with the VLBI core in the VLBI 
images (see Paper I). Note that even sophisticated homo- 
geneous models such as that of H01 do not come close 
to fitting the radio emission, showing that some form of 
inhomogeneous jet component is required. 

While X-ray emission on larger scales in blazar jets 
may be caused by inverse-Compton scattering of the mi- 
crowave background (e.g., Celotti, Ghisellini, & Chiaberge 
2001), on the parsec scales considered in this paper the 
synchrotron photon energy density is orders of magni- 
tude higher than the microwave background energy den- 
sity (about 10 -4 ergs cm -3 compared to 10 -10 ergs cm -3 , 
using expressions from Celotti et al. 2001). Because of 
this large difference in energy densities, inverse-Compton 
scattering of the microwave background is not considered 
further in this paper. 

2. OBSERVATIONS 

The VLBI data used for this paper are taken from the 
22 and 43 GHz VLBI observations of 3C 279 presented in 
Paper I. These data included 18 epochs spanning the time 
range from 1991 to 1997 (at 22 GHz) and 1995 to 1997 (at 
43 GHz). Earlier epochs in this sequence used the Global 
VLBI Network and Mark II recording, later epochs used 
the NRAO VLBA telescopes and correlator. The reader 
is referred to Paper I for further discussion of the VLBI 
observations, and the VLBI images, model fits, and com- 
ponent identifications. Two mosaics at 22 and 43 GHz 



3 



from that paper are shown in Figures 1 and 2 for ease in 
identifying features. Throughout the rest of this paper, 
C4 refers to the moving component approximately three 
milliarcseconds from the core, C5 refers to the stationary 
component at one milliarcsecond from the core, and C5a, 
C6, C7, C7a, C8, and C9 refer to the moving components 
that sequentially emerge from the core during the course 
of the monitoring. See Paper I for more details. High reso- 
lution VLBI data at low frequencies is needed to constrain 
the optically thick VLBI component spectral indices, so we 
also make use of the 1.6 and 5 GHz space VLBI (VSOP) 
observations of 3C 279 from Piner et al. (2000). 

Multiwavelength spectra of 3C 279 are needed to con- 
strain inputs to both the homogeneous sphere and inho- 
mogeneous jet models. For this purpose, we use the set 
of contemporaneous multiwavelength spectra compiled by 
HOI. Eleven such spectra spanning the years 1991 to 2000 
are presented in that paper. We fit these spectra using four 
broken power laws applied to the different regions of the 
spectrum, to provide the appropriate inputs to the mod- 
els. TableQJlists our power-law fits to the multiwavelength 
spectra of H01, for epochs from 1991 to 1997. We exclude 
HOl's epochs P3a (which has no X-ray data) and P5b (the 
very large flare from 1996 February). Spectral indices fol- 
low the naming convention of Konigl (1981): a s i is the op- 
tically thick synchrotron spectral index below the turnover 
frequency, a S 2 is the synchrotron spectral index above the 
turnover frequency but below the break frequency, a S 3 
is the synchrotron spectral index above the synchrotron 
break frequency, and a C 2 is the inverse-Compton spectral 
index in the X-ray region of the spectrum. The sign con- 
vention for spectral indices in this paper is S oc v +a . Er- 
rors on the fitted spectral indices were calculated from 
the errors in the measured flux densities given in H01. 
When only one X-ray measurement was available at a 
given epoch in H01, the X-ray spectral index was acquired 
from the original X-ray paper referenced by H01. Tabled 
also gives the fitted turnover frequency and flux density 
at the turnover frequency (u t and St), the fitted break fre- 
quency and flux density at the break frequency [y^ and 
Sb), and a reference X-ray measurement specified by v x 
and S x that we attempt to match with the models. Figure 
3 shows the power- law fits to the spectra from Tabled 

3. CALCULATION OF SSC DOPPLER FACTORS 

The VLBI morphology and multiwavelength data can- 
not be adequately explained by either an inhomogeneous 
jet (§ 13. 2f) or a homogeneous sphere (or spheres) model 
(§ [EH alone. In a combined model (§ I3.3fl . we require 
the Doppler factor to be the same for the jet and spheres; 
the spheres can therefore be regarded as approximations 
to dense clumps propagating along with the bulk jet. 

3.1. Homogeneous Sphere Model 

Possibly the two simplest blazar models are the homoge- 
neous sphere model (Gould 1979) and the inhomogeneous 
jet (Blandford and Konigl 1979; Konigl 1981). In both 
models, the source is assumed to be moving relativistically 
at a small angle to the line-of-sight, and the electrons re- 
sponsible for the synchrotron emission have a power-law 
distribution in energy. As implied by their names, the ho- 
mogeneous sphere model assumes that the photon density 
within the spherical source is uniform, while the inhomo- 



geneous jet model assumes that the source is a conical jet 
in which the magnetic field and electron number density 
have power-law dependences on the distance from the apex 
of the jet. 

The homogeneous sphere model, the simpler of the two, 
requires six observables to determine the Doppler factor 
and the angle between the source velocity and the line-of- 
sight. Three of these characterize the synchrotron spec- 
trum: the spectral index of the optically thin spectrum, 
and the flux density and frequency at which the sphere 
first becomes optically thick (the "turnover" flux density 
and frequency). Adding the angular size should in prin- 
ciple determine the entire spectrum of the source, if the 
X-rays are the result of SSC scattering of the synchrotron 
photons. In practice, such a calculation predicts X-ray flux 
densities that are several orders of magnitude larger than 
observed, if the bulk motion of the source is assumed to 
be non-relativistic. This "Compton catastrophe" can be 
avoided if the source has a relativistic bulk velocity, and in 
this case the observed X-ray flux density (the fifth observ- 
able) provides a constraint that determines the Doppler 
factor. Together with the apparent superluminal speed 
(the sixth observable), this determines both the bulk ve- 
locity and the angle to the line-of-sight. 

Figure 4 shows sample homogeneous sphere spectra 
plotted with five-frequency VLBI data from late 1997 and 
early 1998. The data at 15, 22, and 43 GHz are from the 
VLB A observation on 1997 Nov 16 (Paper I), while the 
data at 1.6 and 5 GHz are from the VSOP observation 
on 1998 January 9 (Piner et al. 2000). The three bright- 
est components in the source at that epoch (the core, the 
long-lived bright component C4, and the newly emerged 
components C8/9, blended except at 43 GHz) are plotted 
on this figure. Homogeneous sphere spectra are calculated 
from 

S{v) = S {v t /v a T(v/v t ) 2 *{l - O, (2) 

where v t is the turnover frequency, a is the optically thin 
spectral index, So is the flux density at frequency i>o, and 
r = (v t / v) 2 - 5 ~ a . The sphere spectra have been adjusted 
to go through the three highest frequency observations for 
each component. Inspection of Figure 4 shows that homo- 
geneous spheres cannot fit the spectra of either the core or 
the jet components. This is because a homogeneous sphere 
has an optically thick spectral index of 2.5, while the actual 
components have optically thick indices much less than 
this (~ 1.0). Evidently there are inhomogeneities in the 
"components" that broaden the self-absorption turnover 
and cause the spectral index below the turnover to be less 
than 2.5. The high-resolution space VLBI data at low- 
frequencies was crucial to obtaining this result; ground- 
based images at low frequencies blend the components so 
a spectral dissection like that in Figure 4 cannot be ob- 
tained. The primary advantage of the homogeneous sphere 
model is its simplicity; it requires a minimum number of 
observables to constrain the Doppler factor. For this rea- 
son, we proceed with the homogeneous sphere calculation 
in this section, with the expectation that corrections due 
to inhomogeneities (at least for the jet components) will 
be small. An alternative geometry for the VLBI core (a 
conical jet) is tried in the next subsection. An alternative 
geometry for jet components (oblique shocks) is considered 
by, e.g., Aller et al. (2001). 
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Fig. 1.— Time-series mosaic of a selection of 22 GHz VLBI images of 3C 279. Epochs 1991 Jun 24, 1993 Feb 17, 1994 Sep 21, 1996 May 13, 
and 1997 Nov 16 are shown. The images have been restored with a circular 0.2 mas beam without residuals and rotated 25° counterclockwise. 
The lowest contour is 25 mjy beam ; subsequent contours are a factor of two higher than the previous contour. The solid line indicates the 
position of the presumed stationary core. The dotted lines represent the best fits to the model-fit Gaussian positions vs. time. From right to 
left these lines represent C4, C5, C5a, C6, and C7. Some lines have been extended before and after model-fit detections to show speculative 
zero-separation epochs and later positions. 



5 



995.0 



1995.5 



996.0 



to 

CD 



996.5 



1997.0 



997.5 



1998.0 




i | 3C 279 

o parsecs 43 GHz 





o 




o 





oi 






o 



■3 



-4 



mas 



Fig. 2.— Time-series mosaic of a selection of 43 GHz VLBI images of 3C 279. Epochs 1995 Mar 19, 1996 Jan 7, 1996 Jun 9, 1997 Mar 29, 
and 1997 Nov 16 are shown. The images have been restored with a circular 0.15 mas beam without residuals and rotated 25° counterclockwise. 
The lowest contour is 25 mjy beam ; subsequent contours are a factor of two higher than the previous contour. The solid line indicates the 
position of the presumed stationary core. The dotted lines represent the best fits to the model-fit Gaussian positions vs. time. From right to 
left these lines represent C4, C5, C6, C7, C7a, C8, and C9. Some lines have been extended before model-fit detections to show speculative 
zero-separation epochs. 



Table 1 

Spectral Indices and Breakpoints Derived from Multiwavelength Campaign Data, and Reference X-Ray Measurements 



Date Range 


Epocht 


a* 


a 1 


a 1 


a 1 


(Hz) 


(Jy) 


(Hz) 


(Jy) 


(kcV) 


(/"Jy) 


1991 Jun 15 - 


1991 Jun 28 


PI 


0.22 ±0.01 


-0.57 ±0.05 


-1.29 ± 0.06 


-0.67 ±0.03 


6.8 ±3;| xlO 10 


18.6±»;| 


1.1 ±g;§ xlO 13 


i n+ os 


10.0 


0.90 


1992 Dec 22 - 


1993 Jan 12 


P2 


0.24 ±0.02 


-0.81 ±0.03 


-1.66 ± 0.03 


-0.91 ±0.07 


6.5 ±8;| xlO 10 


20.3±°;§ 


1.8 =fcg;| xlO 13 


u - z±3= 0.05 


1.0 


1.00 


1993 Dec 13 - 


1994 Jan 3 


P3b 


0.33 ±0.02 


-0.59 ±0.07 


-1.28 ± 0.03 


-0.75 ±0.17 


5.9±° : | xlO 10 


27.3±J;g 


2.1 ±8;| xio 12 


3 4± 13 


2.3 


0.62 


1994 Nov 29 — 


1995 Jan 10 


P4 


0.22 ±0.01 


-0.63 ± 0.05 


-1.61 ± 0.17 


-0.65 ±0.11 


4.2 ±° l ; 2 xlO 10 


21.2±g;| 


5.2 ±5;° XlO 12 


i n _|_0.8 


2.4 


0.49 


1996 Jan 16 - 


1996 Jan 30 


P5a 


0.33 ±0.03 


-0.50 ±0.12 


-1.72 ± 0.10 


-0.66 ±0.02 


5.7±\-f xlO 10 


26.5± 2 ; 2 


i.i ±8:1 xio 13 


1 Q+1.8 


10.0 


0.33 


1996 Dec 10 - 


1997 Jan 28 


P6a 


0.33 ±0.03 


-0.73 ±0.08 


-1.53 ± 0.07 


-0.67 ±0.02 


4.2 ±8:1 xl ° 10 


28.9±J;g 


i.4 ±J:e xl ° 13 


n 4 Q_l_0.44 


9.5 


0.23 


1997 Jun 17 - 


1997 Jun 24 


P6b 


0.39 ±0.01 


-0.73 ±0.06 


-1.99 ± 0.20 


-0.70 ±0.10 


4.3 ±° 1 ; 2 . xlO 10 


33.8±g;| 


8.3 i^g xlO 13 


14+011 


3.0 


0.48 



f Epoch identification in HOI. 

% a s i is the optically thick synchrotron spectral index below the turnover frequency, a s 2 is the synchrotron spectral index above the turnover frequency but below the break frequency, 
a s 'i is the synchrotron spectral index above the synchrotron break frequency, and a c 2 is the inverse-Compton spectral index in the X-ray region of the spectrum. 

* ut and St are the fitted synchrotron turnover frequency and flux density at the turnover frequency. 

o ut and St are the fitted synchrotron break frequency and flux density at the break frequency. 

o u x and S x are the frequency and flux density of the reference X-ray measurement. 
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3. — Fits of the multiwavelength spectra from Table^with four power laws. The flux density scale in Janskys is correct for the top 
(PI) spectrum, after that each successive spectrum is offset by two orders of magnitude in flux density. The error bars are smaller than the 
plotting symbols. The parameters of the fits are given in Table HI and the epoch designations from TableHlare listed to the right of the figure. 
The spectral indices and breakpoints are labeled for the PI spectrum. 
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Fig. 4. — Five-frequency VLBI spectra of 3C 279 core and jet components from observations between 1997 November and 1998 January. 
Three bright components are shown: the VLBI core (triangles), the blended component C8/9 (diamonds), and C4 (asterisks). The single-dish 
spectrum from Michigan and Metsahovi monitoring is also shown (plus signs). Three sample homogeneous sphere spectra that match the 
three highest frequencies for each component are also shown. 



The Doppler factor of a homogeneous sphere that pro- 
duces SSC X-rays of flux density S x Jy at frequency v x (in 
keV) is (Ghisellini et al. 1993): 



S = f(a)S t 



\n{ubfv t )v% 



6— 4a 5— 3a 



l/(4-2ce) 



(1+z), (3) 



where /(a) « —0.08a + 0.14, St is the flux density (in 
Jy) at the turnover frequency v t (in GHz) obtained by ex- 
trapolating the straight-line optically thin slope (Marscher 
1987), vjj is the synchrotron break frequency in GHz, ^ 
is the angular diameter of the VLBI component (in mas), 
and a is the optically thin spectral index. Equation J3J 
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finds the Doppler factor for which the given spherical com- 
ponent produces all the X-ray flux density from the source. 
Since the other components must also contribute X-rays 
at some level, the Doppler factor found by this method is 
a lower limit for the given component. To put this method 
into practice to find an actual Doppler factor for 3C 279, 
we must assume the source has a constant Doppler fac- 
tor along its jet. We then apply equation © to find the 
lower limit to the Doppler factor for each component at 
a given epoch. The actual Doppler factor at that epoch 
must then be somewhat larger than the highest lower limit 
found from equation J2J. We then calculate the Doppler 
factor that will produce the observed X-ray flux density 
from the sum of the X-rays from all components, and what 
percentage of the observed X-ray flux density is produced 
by each component. 

In Table |21 we show the inputs to equation © and the 
resulting Doppler factor lower limits at each 22 GHz VLBI 
epoch, for those components found to produce a significant 
percentage (> 1%) of the X-ray flux density at that epoch. 
The parsec-scale radio structure of 3C 279 consists of the 
compact core, the bright long-lived component C4 located 
about 3 mas from the core at these epochs, and a series 
of short-lived inner-jet components that have all faded by 
the time they reached about 1 mas from the core. The 
inner-jet components have similar spectra, typified by the 
spectra of the blended component C8/9 in Figure 4. We 
use C8/9 as a guide for all of the inner-jet components, 
and take them to have v t = 22 GHz and a = —0.6 (see 
Figure 4). For component C4, we take i> t = 7 GHz and 
a = —0.5, based on power-law fits to the optically thick 
and thin portions of the spectrum in Figure 4. For the core 
turnover frequency and spectral index, we use the turnover 
frequency and spectral index (0^2) of the multiwavelength 
spectrum from the HOI epoch closest in time to the VLBI 
epoch (see Table QJ, since the core dominates the high- 
radio-frequency spectrum. St for the jet components is 
calculated from the observed flux density at the turnover 
frequency (for the inner-jet components), or extrapolated 
to 7 GHz (for C4), with the opacity correction described 
by Marscher (1987) applied (St = S b s e Tm ), where r m is 
tabulated by Marscher (1987). The core St is taken to be 
the greater of St from the closest multiwavelength spec- 
trum with extrapolations of jet component flux densities 
subtracted, or the observed 22 GHz (or 43 GHz if avail- 
able) VLBI core flux density with the opacity correction 
applied. £<j is taken to be the VLBI model-fit size (or size 
upper limit if the fitted size is zero), multiplied by 1.8 to 
convert from the model fit Gaussian FWHM to the diame- 
ter of an optically thin sphere (Pearson 1995). Uf,, v x , and 
S x are taken from Tabled from the multiwavelength spec- 
trum closest in time to the VLBI epoch. Note equation ® 
is very insensitive to these three values (e.g., 5 oc S x 2 for 
a — —0.5), X-ray variability between the X-ray and VLBI 
observation should not have a large affect (e.g., if the X- 
ray flux density increases by a factor of 10, the calculated 
Doppler factor increases by a factor of 1.6). 

In Table we calculate the Doppler factor for each 
epoch, assuming that all components at a given epoch 
have the same Doppler factor, and that the total X-ray 
emission from all components must equal the observed X- 
ray emission. This Doppler factor will be at least slightly 
higher than the highest lower limit found for that epoch 



in Table |2 since each lower limit in Table [5] assumed the 
component under consideration was the sole source of the 
X-ray emission. Errors in the turnover frequency and an- 
gular size in equation cause large errors in the calcu- 
lated Doppler factor (8 oc v\ z and £^' 6 f° r a — ~ 0.5). 
Because our measurements of these quantities are only ac- 
curate to about 25%, we estimate our calculated Doppler 
factor lower limits in Tables [21 and [31 to be correct only to 
within a factor of two. The Doppler factor values listed 
in Table 13 confirm this: the average Doppler factor is 23, 
and the measured values show about a factor of two scat- 
ter around this value, with most falling between 10 and 
40. There is some indication that the Doppler factor in- 
creases with time from 1991 to 1997, but given the errors in 
the calculated Doppler factor, this may not be significant. 
The calculated relative X-ray brightness of components 
in turn depends on the Doppler factor from Table [5] to a 
high power (S x oc S 5 for a — —0.5), so the relative X-ray 
brightnesses of components can be in error by as much as 
a factor of 30. This means that any of the components 
listed as producing more than several percent of the X-ray 
emission at a given epoch in Tablets a candidate for pro- 
ducing most of the X-rays. At eight of the epochs listed in 
Table the calculations show that the core produces at 
least two orders of magnitude more X-rays than the next 
brightest component, and the calculations show the core to 
be the leading X-ray producer at 15 of the 18 epochs. De- 
spite the large errors in the calculations, we can conclude 
that in the homogeneous sphere model, the VLBI core is 
the dominant source of the X-ray emission from late 1992 
through 1997. It appears that C4 dominated the X-ray 
emission for a time prior to 1992 (due to its small size at 
that time), and that there may also be contributions at 
about the 10% level from whatever inner-jet component 
has just emerged from the core (this changes from C6 to 
C9 over the course of the observations) . 

3.2. Inhomogeneous Conical Jet Model for Core 

In this subsection we consider an alternate geometry for 
the VLBI core, the inhomogeneous jet model of Konigl 
(1981). The jet is represented by a cone with opening 
half-angle cf>, and the axis of the jet makes an angle 9 
with the line-of-sight (6 > 4>). The bulk Lorentz factor 
of the jet is 7, the electron Lorentz factor is j e . Electron 
Lorentz factors lie between the limits 7 e ; and 7 eu , where 
j e i is set to 100 for all models. The magnetic field and 
electron number density (or electrons plus positrons) are 

determined by B = B\r~ m and n e (^ e ) = K e r~ n ji 2ae 1 \ 
where r is the distance in parsecs from the apex of the 
jet, a e — a and K e = Kl for j e < 7 e h, and a e = a — 0.5 
and K e — Klj e b for j e > 7 e fc, where 7 e & is a function 
of r (j e i < 7 e h < 7e«). The location of the break in the 
power-law is estimated by equating the jet travel time to 
a distance r with the synchrotron cooling time; see equa- 
tion (4) of Konigl (1981) and equation (21) of Blandford 
& Konigl (1979). The VLBI 'core' emission is then due to 
the integrated spectrum of the unresolved conical jet, and 
the position and size of the core change with frequency 
in a predictable way. At low frequencies, the dominant 
emission region follows the local turnover frequency and 
moves in with increasing frequency (r oc i/ -1 /*^ where 
k m = ((3 — 2a)r7i+2n — 2)/(5 — 2a)), and the spectral index 
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Table 2 

Lower Limits to the Doppler Factor Computed Using the Homogeneous Sphere Model 











s t 






u 




VLBI Epoch 


HOI Epoch 


Comp. 


(Jy) 


(GHz) 


a 


(mas) 


5 


1991 Jun 24 


PI (1991 Jun 15 - 


1991 Jun 28) 


C4 


4.1 


7.0 


-0.5 


0.25 


15.7 


1992 Jun 14 


P2 (1992 Dec 22 - 


1993 Jan 12) 


Core 


20.5 


65.0 


-0.81 


0.20 


5.9 








C4 


4.4 


7.0 


-0.5 


0.49 


7.3 


1992 Nov 10 


P2 (1992 Dec 22 - 


1993 Jan 12) 


Core 


18.6 


65.0 


-0.81 


0.14 


10.0 




C6 


4.4 


22.2 


-0.6 


0.25 


4.2 








C4 


6.2 


7.0 


-0.5 


0.83 


4.4 


1993 Feb 17 


P2 (1992 Dec 22 - 


1993 Jan 12) 


Core 


22.7 


65.0 


-0.81 


0.16 


9.1 








C4 


5.2 


7.0 


-0.5 


0.79 


4.0 


1993 Nov 8 


P3b (1993 Dec 13 - 


- 1994 Jan 3) 


Core 


24.7 


58.7 


-0.59 


0.16 


11.8 


1994 Mar 2 


P3b (1993 Dec 13 - 


- 1994 Jan 3) 


Core 


23.5 


58.7 


-0.59 


0.14 


13.8 








C7 


6.7 


22.2 


-0.6 


0.20 


8.7 


1 AA /IT 1 A 

1994 Jun 12 


P4 (1994 Nov 29 - 


1995 Jan 10) 


Core 


28.0 


42.4 


-0.63 


A 1 O 

0.18 


1 D T 

18.7 


1994 Sep 21 


P4 (1994 Nov 29 - 


1995 Jan 10) 


Core 


26.2 


42.4 


-0.63 


0.18 


17.5 


1995 Jan 4 


P4 (1994 Nov 29 - 


1995 Jan 10) 


Core 


21.8 


42.4 


-0.63 


0.13 


25.9 


1995 Feb 25 


P4 (1994 Nov 29 - 


1995 Jan 10) 


Core 


17.6 


42.4 


-0.63 


0.16 


14.0 








C7 


3.1 


22.2 


-0.6 


0.11 


11.8 


1995 Mar 19 


P4 (1994 Nov 29 - 


1995 Jan 10) 


Core 


22.7 


42.4 


-0.63 


0.14 


21.8 


1 a a /"> t >"7 

1996 Jan 7 


P5a (1996 Jan 16 - 


- 1996 Jan 30) 


Core 


23.2 


56.6 


-0.5 


0.11 


26.1 




C4 


4.2 


7.0 


-0.5 


0.40 


9.6 


1996 May 13 


P5a (1996 Jan 16 - 


- 1996 Jan 30) 


Core 


24.3 


56.6 


-0.5 


0.11 


27.2 




C7a/C8 


8.7 


22.2 


-0.6 


0.20 


11.4 


1996 Jun 9 


P5a (1996 Jan 16 


- 1996 Jan 30) 


Core 


25.3 


56.6 


-0.5 


0.07 


54.4 


1997 Jan 15 


P6a (1996 Dec 10 - 


- 1997 Jan 28) 


Core 


26.2 


42.4 


-0.73 


0.12 


30.1 




PR 


14.0 


00 

zz.z 


-U.D 


n oo 
u.zz 


17 








C4 


6.6 


7.0 


-0.5 


0.31 


24.6 


1997 Mar 29 


P6a (1996 Dec 10 - 


- 1997 Jan 28) 


Core 


20.7 


42.4 


-0.73 


0.10 


31.0 






C8/9 


16.1 


22.2 


-0.6 


0.18 


26.8 








C4 


4.6 


7.0 


-0.5 


0.20 


34.5 


1997 Jul 16 


P6b (1997 Jun 17 


- 1997 Jun 24) 


Core 


24.8 


42.8 


-0.73 


0.13 


26.8 




C8/9 


16.9 


22.2 


-0.6 


0.31 


12.4 








C4 


7.8 


7.0 


-0.5 


0.43 


17.0 


1997 Nov 16 


P6b (1997 Jun 17 


- 1997 Jun 24) 


Core 


41.4 


42.8 


-0.73 


0.14 


38.1 




C4 


8.9 


7.0 


-0.5 


0.40 


22.3 



of the integrated spectrum is a s i- At radius ru and fre- 
quency v t , the local break frequency becomes less than the 
turnover frequency, and the dominant emission region be- 
gins to follow the local break frequency and moves out with 
increasing frequency (r oc v 1 / kb J where kb = 3m — 2), until 



7 eu at frequency and radius r = (yb/^t) 



i/feb 



which we set to be r u , the upper radius of the Konigl jet. 
We use this set of conditions to determine "f eu in our mod- 
els. The spectral index of the integrated spectrum in this 
region is a s2 . Above i>b, the dominant emission region 
moves in with frequency, following the increasing mag- 
netic field (r oc v^ 1 / m ), and the integrated synchrotron 
spectrum falls with index a s z- We extend the calculated 
emission from the upper synchrotron branch interior to 
Tm, and that from the lower synchrotron branch exterior 
to r u , rather than truncating the spectrum at these radii 
as Konigl does. The reader is referred to Hutter & Muf- 
son (1986) for an illuminating diagram of this model in 
the radius vs. frequency plane. 



3.2.1. Calculation of a, m, and n 

The parameters a, to, and n completely determine the 
synchrotron and SSC spectral indices of the integrated 
spectrum by the following equations (Konigl 1981): 



"si 



-15 + 4a + 5m — 4am + 5n 
3to - 2am + 2n - 2 



a.s2 



max 



4am-2a-m-n+3 
3m-2 

a - 0.5 



a s3 = -(to + 2 - n)/m 



a c i = mm 



-4am-2an+6a+5m+10n-20 
3m — 2am+2n- 2 
1 



a c2 = a - 



(1 - a)m + 2n - 4 
7m - 4 



(4) 

(5) 
(6) 

(7) 
(8) 
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Table 3 

X-Ray Production in the Homogeneous Sphere Model 

Fraction of Total X-Rays by Component 

C6 C7 C7a/C8t C8 08/9* 



0.01 



Epoch 


5 


Core 




C4 


1991 


Jun 24 


16 






1 


00 


1992 


Jun 14 


8 





22 





78 


1992 


Nov 10 


10 





97 





02 


1993 


Feb 17 


9 





98 





02 


1993 


Nov 8 


12 


1 


00 






1994 


Mar 2 


14 





92 






1994 


Jun 12 


19 


1 


00 






1994 


Sep 21 


18 


1 


00 






1995 


Jan 4 


26 


1 


00 






1995 


Feb 25 


15 





71 






1995 


Mar 19 


22 


1 


00 






1996 


Jan 7 


26 





99 





01 


1996 


May 13 


27 





99 






1996 


Jun 9 


54 


1 


00 






1997 


Jan 15 


32 





69 





26 


1997 


Mar 29 


39 





30 





55 


1997 


Jul 16 


27 





89 





09 


1997 


Nov 16 


39 





94 





06 




0.15 
0.02 



t Refers to a blended component. 

Table 4 

Fitted values of Conical- Jet Spectral Parameters 



Epoch 


a 


m 


n 


a s i 


Os2 


a S 3 


0L C 2 


h 


PI 


-0.4 


1.6 


1.5 


0.21 


-0.66 


-1.31 


-0.57 


1.2 


P2 


-0.6 


2.0 


0.8 


0.18 


-0.85 


-1.60 


-0.68 


1.3 


P3b 


-0.3 


1.8 


1.5 


0.33 


-0.55 


-1.28 


-0.46 


1.3 


P4 


-0.4 


1.7 


1.4 


0.22 


-0.65 


-1.35 


-0.55 


1.3 


P5a 


-0.5 


2.0 


1.1 


0.30 


-0.78 


-1.45 


-0.62 


1.4 


P6a 


-0.5 


1.9 


1.3 


0.34 


-0.81 


-1.37 


-0.66 


1.4 


P6b 


-0.5 


1.9 


1.4 


0.39 


-0.84 


-1.32 


-0.68 


1.4 



The SSC index a c \ (associated with the optically thick 
portion of the synchrotron spectrum) is not observable in 
the multiwavelength spectrum because the SSC emission 
in this frequency range lies well below the synchrotron 
emission. The SSC index a C 2 (associated with the op- 
tically thin portion of the synchrotron spectrum) is the 
SSC index that is actually observed in the X-ray portion 
of the multiwavelength spectrum. 

Since the Konigl model parameters are so closely tied to 
the multiwavelength spectrum, we make one calculation of 
S for each multiwavelength epoch listed in Tabled rather 
than one for each VLBI epoch as in § 13.11 We fit values 
of a, to, and n at each epoch by finding the combination 
that most closely reproduces the four observable spectral 
indices, given in Table [T]for each multiwavelength epoch. 
We take into account the errors on the observed spectral 
indices given in Tabled an d find the values of a, to, and 
n that minimize x 2 a t each epoch. 



Fitted values of a, m, and n are given in Table 0] along 
with the spectral indices calculated from these best-fit val- 
ues, for comparison with Tabled In some cases, only mod- 
est agreement could be made with some of the observed 
indices. Note that the observed a s \ is an integrated value 
for the source that includes some steep-spectrum jet emis- 
sion. The observed a s \ is therefore a only limiting value 
for the core a s \, and later we modify the fits so the ra- 
dio emission from the Konigl model lies slightly below the 
total observed radio emission. The magnetic field index 
to is restricted to lie between 1 (for a purely transverse 
field) and 2 (for a purely longitudinal field). Fitted values 
for a lie between —0.3 and —0.6, those for to between 1.6 
and 2.0, and those for n between 0.8 and 1.5. In Table 01 
we also give Konigl's parameter k m that controls the fre- 
quency dependent size of the core, which is proportional to 
j,-l/fc"\ This provides an independent check on the model 
fits, and using the mean core sizes at 22 and 43 GHz from 
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Paper I we calculate a k m of 1.2, in good agreement with 
the mean k m of 1.3 from Table 

3.2.2. Calculation of Jet Doppler Factor 

The fitted values of a, to, and n from Table 0] along 
with other observable parameters, were used as inputs to 
a Konigl model calculation at each epoch, implemented 
in the Mathcad software package. The other input pa- 
rameters were: the synchrotron turnover frequency and 
flux density and the synchrotron break frequency from Ta- 
ble n the observed superluminal speed, and the projected 
distance of the VLBI core from the apex of the jet, r pro j. 
The angle of the jet to the line-of-sight was then varied 
until the predicted X-ray flux density from the Konigl jet, 
including both synchrotron and SSC flux density, matched 
the observed X-ray flux density S x at frequency v x from 
Tabled Once the angle to the line-of-sight is known, the 
observed superluminal speed then determines the intrinsic 
jet speed and Doppler factor, and the opening half-angle 
of the jet cf> can be calculated from equations in Konigl 
(1981). 

Input parameters not taken from Table ^ were deter- 
mined as follows. The superluminal speed was taken to 
be the measured speed of the inner-jet components C6, 
C7, C7a, and C8 (the components closest to the core dur- 
ing the relevant epochs) from Paper I. Measured speeds of 
these four components are consistent with a constant aver- 
age inner-jet speed of 5.2c (see Table 4 of Paper I), which 
was the apparent speed adopted for these model fits. The 
projected distance of the VLBI core from the apex of the 
jet, r pro j, is not a directly observable quantity, but does 
effect two observable quantities: the measured core radius 
£ r rj r pro j sine/)/ sin 8, and the apparent jet opening angle 
<t>app ~ <fi/ sin#. In practice, r pro j was varied until the best 
fit was obtained to the observed values of these two quan- 
tities: £ r = 0.059 ± 0.018 mas (using the mean FWHM of 
the Gaussian model fits from Paper I, and multiplying by 
1.6 to convert to the diameter of an optically thick sphere 
[Pearson 1995]), and 4> a pp = 10.5 ± 5.3° (using the mean 
opening angle given by component C4 from Paper I, with 
C4's Gaussian FWHM converted to the diameter of an 
optically thin sphere for consistency). 

Results from this model fitting, including the adopted 
value of r pro j and the calculated values of 4> a pp, £r, 0, 5, 
and 7, are given in Table 03 Also given in Table are 
tm (the smallest radius from which synchrotron emission 
with index a is observed), r u (in our implementation, the 
largest radius from which synchrotron emission with index 
a — 0.5 is observed), the magnetic field and relativistic- 
particle number density at tm, and the ratio of relativistic 
particle to magnetic energy density at tm and r u . A sketch 
of a sample geometry from Table [S] is shown in Figure 5. 
The model synchrotron and SSC spectra are shown in Fig- 
ure 6, where the sums of these spectra are compared with 
the observed multiwavelength spectra. 

In reality, some of the radio emission from 3C 279 is due 
to the parsec-scale jet that is resolved by VLBI observa- 
tions, and not to the partially resolved core. To take this 
into account, we constructed an alternate set of model fits 
with v t increased by a factor of 1.5 and St decreased by 
a factor of 0.75 (relative to their values in Table 0). This 
causes the predicted core radio flux density to lie slightly 
below the observed single-dish radio flux densities, allow- 



ing for some emission from the extended jet (Piner et al. 
2000; de Pater & Perley 1983), which does not contribute 
significantly to the integrated spectrum above the turnover 
frequency. In these models the jet contributes ~ 3 Jy at 
5 GHz, and these model fits are indicated by an asterisk 
next to the epoch name in Table 03 This set of models is 
used in all subsequent calculations. 

The agreement between the calculated and observed 
spectra is reasonably good, although there are discrepan- 
cies (particularly with the optical flux densities, e.g., 1994 
Nov 29 — 1995 Jan 10, period P4 in H01) of about a factor 
of 2. There are several reasons for this: the optical flux 
density of 3C 279 is known to vary by as much as 60% 
on timescales of a day (Balonek & Kartaltepe 2002) so 
variability may be a factor, in some cases one of the four 
observed spectral indices was poorly matched by the three 
free parameters a, to, and n, and at some epochs there are 
apparently significant synchrotron X-rays, so the observed 
a C 2 is not an accurate measurement of the Compton spec- 
tral index at those epochs. The spectra in Figure 6 do not 
extend up to the EGRET energy range, and we do not 
attempt to model the high-energy 7-ray emission with the 
Konigl model. Although the Konigl model can produce 
7-rays in the EGRET energy range, there is considerable 
evidence that the GeV emission is external Compton and 
not SSC (Kubo et al. 1998; H01; Ballo et al. 2002), so we 
do not apply the EGRET data as a constraint on our SSC 
model. 

The values of the Doppler factor obtained by applying 
a conical geometry to the core are about a factor of two 
lower than the Doppler factors obtained using a spheri- 
cal geometry (§ I3.1fl . showing that the assumed geometry 
can have a moderate influence on the SSC Doppler fac- 
tor. The average Doppler factor from the fits in Table 
that allow for emission from an extended jet is 9 with a 
scatter of ±3. When inputs to the model are all varied by 
10% (such that the changes all act together to increase or 
decrease the Doppler factor), the fitted 5 can be made to 
vary by about 30%. Thus variations in the Doppler factor 
from epoch to epoch in Table [S] could be due to the ob- 
servational uncertainties in the input parameters, rather 
than real variations in the Doppler factor. 

3.3. Combined Conical- Jet/ Sphere Models 

We are now in a position to make a model for the source 
that combines the conical VLBI core with the homoge- 
neous sphere VLBI jet components (§ I3.1|l . We use the 
homogeneous sphere values from Table El and the conical- 
jet core values from Table using the multiwavelength 
epoch closest in time to each VLBI epoch. As in S I3.1I each 
calculated Doppler factor (either conical jet or sphere) pro- 
vides a lower limit to the Doppler factor at that epoch 
(assuming a constant Doppler factor along the jet), since 
each calculation assumes the component under considera- 
tion produces all of the X-rays. We solve for one Doppler 
factor at each epoch as follows. We lower the X-ray flux 
density attributed to the conical-jet core in the Konigl-jet 
fitting program, while adjusting r pro j to maintain the ob- 
served values of 4> a pp and £ r . This has the effect of raising 
the calculated Doppler factor. We use this Doppler factor 
to calculate the X-ray flux density that would be observed 
from the jet components in Table |2 and add these flux 
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Table 5 

Results from Conical- Jet Model Fits 











n 
V 






TM 


r u 


nt 


n e T 


TT ITT { \ t 


t t It J / \ t 

Ue/U B {r u y 


Epoch* 


(mas) 


(deg) 


(mas) 


(deg) 


7 


c 



(PC) 


(PC) 


(G) 


(cm s ) 


PI 


0.27 


11.6 


0.055 


14.0 


5.5 


4.0 


2.6 


15.8 


0.23 


42 


5.1 


338 


PI* 


0.18 


14.4 


0.045 


11.8 


5.3 


4.9 


1.5 


7.7 


0.30 


73 


4.8 


234 


P2 


0.28 


11.8 


0.058 


6.5 


5.8 


8.1 


6.3 


25.7 


0.17 


5.5 


0.9 


168 


P2* 


0.19 


13.6 


0.045 


4.9 


6.3 


9.8 


4.1 


15.2 


0.22 


9.7 


0.9 


122 


JroD 


n 01 




U.U45 


Q Q 
O.O 


7 Q 
1 .0 


1 O Q 

Iz.o 




ZO.O 


fl 17 


1 Q 
1 .0 


n q 
U.o 


D 


P3b* 


0.16 


11.7 


0.033 


2.6 


8.1 


14.3 


7.4 


18.6 


0.21 


3.3 


0.4 


7 


P4 


0.31 


10.9 


0.059 


10.0 


5.3 


5.7 


6.2 


29.1 


0.16 


4.2 


0.8 


142 


P4* 


0.22 


11.4 


0.044 


7.8 


5.5 


7.0 


4.0 


16.7 


0.21 


8.5 


1.1 


51 


P5a 


0.36 


9.4 


0.059 


7.8 


5.5 


7.1 


7.8 


29.0 


0.14 


7.6 


2.0 


223 


P5a* 


0.23 


12.9 


0.051 


6.3 


5.8 


8.3 


4.5 


15.1 


0.19 


12 


1.7 


126 


P6a 


0.50 


7.6 


0.066 


10.0 


5.3 


5.7 


10.4 


49.4 


0.12 


6.8 


2.7 


343 


P6a* 


0.31 


10.5 


0.057 


8.2 


5.5 


6.8 


5.8 


24.8 


0.16 


11 


2.2 


205 


P6b 


0.45 


8.4 


0.066 


5.1 


6.2 


9.5 


18.4 


143 


0.10 


1.5 


0.7 


323 


P6b* 


0.28 


10.9 


0.053 


3.5 


7.2 


12.1 


12.6 


87.2 


0.14 


2.0 


0.5 


173 



*Models indicated by an asterisk have vt increased by a factor of 1.5 and St decreased by a factor of 0.75 relative to their values in Table HI 
This allows for some emission from the extended jet, which causes the core radio flux density to lie slightly below the observed single-dish radio 
flux density. 

t Magnetic field and total electron number density are evaluated at radius tm ■ 
* Ratio of electron to magnetic energy density at tm and r u . 



densities to that attributed to the VLBI core. We con- 
tinue this process until we find a total X-ray flux density 
that equals the observed X-ray flux density. Results of this 
process are given in Table El The parameter space of the 
Konigl jet model is complicated, and there is not a simple 
relation between the observed X-ray flux density and 8 as 
for the homogeneous sphere. In most regions of parame- 
ter space 8 is a slowly varying function of the X-ray flux 
density as for a sphere, but there are regions where 8 can 
change dramatically in response to a small change in the 
X-ray flux density (e.g., the last two entries in Table 
where, as the solution for 9 begins to approach zero, the 
Lorentz factor increases to maintain the observed superlu- 
minal speed, producing the high Doppler factors in these 
models). 

We reiterate that the predicted X-ray emission depends 
strongly on the observed quantities, so that any of the 
components listed as producing more than several percent 
of the X-ray emission at a given epoch in Table is a can- 
didate for producing most of the X-rays. Because of this, 
it is better to consider averages over many epochs than the 
entry in Table [5] from a single epoch. On average, mod- 
eling the core as a conical jet rather than a homogeneous 
sphere reduces its contribution to the X-ray flux density. 
When the core is modeled as a conical jet, then on average 
the core produces about half of the X-rays, with the other 
half being split about evenly between C4 and the brightest 
inner-jet component. 

4. DISCUSSION 
4.1. Comparison to Other Estimates of 8 for 3C 279 



A strong lower limit to the Doppler factor can be ob- 
tained by enforcing the condition that the emitting region 
should be transparent to 7-rays. Inferring a source size 
from the 7-ray variability timescale during the large flare 
in early 1996, Wehrle et al. (1998) find 8 > 6.3 for 1 GeV 
photons and 8 > 8.5 for 10 GeV photons. This agrees very 
well with the lower limit to the core 8 found at this epoch 
in this paper, 8 > 8.3 (Table [SJ, and with the overall 8 
found at this epoch, 8=11 (Table|SJ). 

An independent method for measuring 8 is to compare 
the radio core brightness temperatures measured from 
VLBI maps and from radio light curves (Lahteenmaki et 
al. 1999). Since these depend on the intrinsic brightness 
temperature multiplied by 8 raised to different powers, a 
measured value of 8 can be extracted from these two ob- 
servables. We applied this method to 3C 279 in Paper I 
and found <5=7.4 at epoch 1995.2. Again, this agrees very 
well with the lower limit to the core 8 found at this epoch 
in this paper, 8 > 7.0 (Table [SJ, and with the overall 8 
found at this epoch, 8=9 (Table EJ. The similar Doppler 
factors found from these two independent methods dif- 
fer from the much higher Doppler factor (8 ~ 100) found 
when equipartition between magnetic and particle energy 
is assumed (Paper I). Other calculations of an equiparti- 
tion Doppler factor for 3C 279 (e.g. Giiijosa & Daly 1996) 
have found a lower equipartition Doppler factor because 
they did not have access to the high-resolution (and high 
brightness temperature sensitivity) VLBI data presented 
in Paper I. 

Models using various emission processes and geometries 
to explain the multiwavelength spectra of 3C 279, includ- 
ing the 7-ray emission, have constrained 8 as part of their 
model fitting. Some examples are as follows: Maraschi, 



13 



pro: 




( r (VLB core size at 22 GHz) 



ine of sight 



Fig. 5. — Geometry of the conical-jet model, with numerical values taken from the first row of Tablc|S](thc fit to the 1991 Jun 15 — 1991 
Jun 28 (PI) epoch). Note that <f> is the opening half-angle of the jet, £ r is the jet cross-sectional radius, r is the linear distance along the jet, 
and r pro j is the projected linear distance along the jet. 



Ghiscllini, & Celotti (1992) used an accelerating parabolic 
jet where 5 varied from 10 to 18. Ghiscllini & Madau 
(1996) assumed 5=14 for application of their "mirror" 
model to the high-energy emission of 3C 279. HOI and 
Ballo et al. (2002) find Doppler factors ranging from 8 to 
23 and from 12 to 19 respectively for application of their 
models to the various multiwavelength spectra presented 
in those papers. In these cases the angle to the line-of- 
sight and / or the Lorentz factor are assumed input quanti- 
ties, so the derived Doppler factors are merely consistent 
within the framework of the particular model, not actual 
Doppler factor measurements. In addition, these models 
do not apply constraints from the VLBI observations as 
we do in this paper, so they predict features in the VLBI 
maps that are not observed. The models mentioned in this 
paragraph predict apparent speeds in the inner jet rang- 
ing from 1 to 24c, whereas the observations show apparent 
speeds in the inner jet to be about 5c (Paper I), which is 
enforced by all of the models in Table [SJ 

Note that our reliance on the apparent superluminal 
speed assumes that the pattern speed observed in the 
VLBI observations (Paper I) is equal to the bulk fluid 
speed. While there is one stationary component (C5) that 
clearly does not move at the bulk fluid speed, all five com- 
ponents observed in the inner jet of 3C 279 (C5a, C6, 
C7, C7a, and C8) during the course of the monitoring de- 
scribed in Paper I moved with approximately the same 
apparent speed of 5c, so we take this value as an indicator 



of the apparent fluid speed in the inner jet. If the appar- 
ent speed is not used as an observable, then the speed and 
orientation of the jet are not tightly constrained in our 
models. For example, allowing the apparent bulk speed 
to range from 1 to 10c allows the following ranges of pa- 
rameters for the PI model in Table 5° < 6 < 14°, 
3 < 7 < 24, and 2 < 6 < 4. 

4.2. Implications of Konigl Jet Model 

4.2.1. Comparison to Sphere Model Results 

In comparing conical and spherical models for the VLBI 
core, we find that a smaller Doppler factor is required to re- 
duce the predicted core X-ray flux density to the measured 
X-ray flux density using the conical geometry. This implies 
that, if a conical geometry is indeed the correct geometry 
for the VLBI core, inverse- Compton calculations such as 
those of Ghisellini et al. (1993) that have assumed a spher- 
ical geometry may have systematically overestimated 5 in 
their samples. 

4.2.2. Implications of a, m, and n 

The values of a, m, and n derived for the Konigl model 
determine the orientation of the magnetic field and the ra- 
tio of relativistic particle to magnetic energy density. The 
value of m can vary between 1 and 2 in the Konigl model, 
with m = 1 corresponding to a purely transverse mag- 
netic field, and m — 2 to a purely longitudinal field. Con- 
servation of particle number in a conical jet requires that 
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Fig. 6. — Sum of the calculated synchrotron and inverse-Compton spectra from the Konigl model fit to each multiwavelength epoch, 
compared with the observed flux densities at each multiwavelength epoch, indicated by asterisks. The spectral sum is indicated by the solid 
line, the synchrotron spectrum by the dotted line, and the inverse-Compton spectrum by the dashed line. 
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Table 6 

X-Ray Production in the Combined Conical- Jet/Sphere Model 



Fraction of Total X-Rays by Component ffl 
Epoch 5 Core C4 C6 06/7* C7 C7a/C8t C8 Cg/gj (deg) 7* 



1991 Jun 24 


16 





09 





91 










2 


8 


9.8 


1992 Jun 14 


10 





83 





17 










4 


4 


6.6 


1992 Nov 10 


10 





97 





02 


0.01 








4 


8 


6.4 


1993 Feb 17 


10 





99 





01 










4 


9 


6.3 


1993 Nov 8 


14 


1 


00 














2 


6 


8.1 


1994 Mar 2 


15 





95 








... 0.05 






2 


3 


8.6 


1994 Jun 12 


8 





77 





15 




... 0.08 






6 


6 


5.7 


1994 Sep 21 


8 





84 





01 




0.15 






6 


9 


5.7 


1995 Jan 4 


8 





72 





07 


0.05 


... 0.16 






6 


3 


5.8 


1995 Feb 25 


13 





36 





01 




... 0.63 






3 


2 


7.4 


1995 Mar 19 


9 





57 





06 




... 0.37 






5 


3 


6.1 


1996 Jan 7 


11 





48 





51 






0.01 




4 


9 


8.1 


1996 May 13 


12 





28 





16 






0.56 




3 


4 


7.2 


1996 Jun 9 


19 





08 





18 






0.74 




1 


6 


10.2 


1997 Jan 15 


26 





08 





77 






... 0.15 




1 


2 


14.1 


1997 Mar 29 


37 





09 





72 








0.19 





6 


19.3 


1997 Jul 16 


27 





87 





11 








0.02 





8 


13.8 


1997 Nov 16 


33 





85 





15 













5 


16.8 



t Refers to a blended component. 

\ These values are discussed in § 14.2.41 



from spectral fits (rather than assumed). Similar particle 
dominance was found by Unwin et al. (1994, 1997) and 
by Hutter and Mufson (1986), who attributed the result 
qualitatively to a conversion of magnetic energy to parti- 
cle energy by magnetohydrodynamic jet acceleration. If 
the Konigl model is a correct description of the VLBI core 
of 3C 279, then this core is not in equipartition, which 
could explain why 3C 279's equipartition Doppler factor 
is so much higher than the Doppler factor measurements 
obtained by other means (see § 14. 1|) . This finding con- 
trasts with that of Lahteenmaki et al. (1999), who find 
that Doppler factors computed by comparing variability 
and VLBI brightness temperatures in general agree with 
equipartition Doppler factors, although analysis of bright- 
ness temperature measurements from the high-resolution 
VLBI data of Paper I has shown that this is not the case 
for 3C 279. 

From the particle energy density, we can calculate the 
energy flux associated with the particles from 

L kln w ^ircr 2 jet j 2 U e (l + k), (9) 

where Tjet is the linear size of the jet cross-sectional radius, 
U e is the relativistic-particle energy density (see Table |5J, 
and k is the ratio of proton to electron energy. See De 
Young (2002) equation (4.107), Bicknell (1994) equation 
(52), and Celotti & Fabian (1993) equation (1) for dis- 
cussions of this equation. Note that U e depends on 7^, 
which is not precisely known (although the observed lack 
of Faraday rotation in most extragalactic radio sources 
implies j e i > 100, Jones & O'Dell 1977), so that this is 
only an order-of-magnitudc calculation. From equation El 



n = 2, our values of n < 2 imply that the total number 
of relativistic particles increases down the jet (e.g., from 
continuous acceleration adding to the total number of rel- 
ativistic electrons). Our fitted values of a vary between 
—0.3 and —0.6, our values of n between 0.8 and 1.5, and 
our values of m between 1.6 and 2.0, which corresponds to 
a predominantly longitudinal magnetic field in the region 
of the jet modeled by the Konigl model, which is about 
0.1 mas in size for 3C 279. In contrast, VLBI polarime- 
try observations at many frequencies from 15 to 86 GHz 
(Leppanen, Zensus, & Diamond 1995; Taylor 1998; Lister, 
Marscher, & Gear 1998; Homan & Wardle 1999; Lister & 
Smith 2000; Attridge 2001) have given magnetic field vec- 
tors oriented perpendicular to the jet in the core region. 
We do not consider these results to be in conflict, since the 
VLBI polarimetry observations show the core to have low 
levels of polarization even at 86 GHz (Attridge 2001), sug- 
gestive of an initially tangled field that becomes ordered 
at shocks outside the Konigl jet region. VLBI polarime- 
try with beams smaller than 0.1 mas would be needed to 
image the magnetic field structure in the region modeled 
in this paper. 

The values of m and n also determine how the ratio 
of relativistic particle to magnetic energy density varies 
along the jet (it goes approximately as r ~ n+2m ), this ratio 
is quoted in Table|S]at radii tm and r u . The jet is close to 
equipartition at Tm, but the degree of particle dominance 
increases down the jet, and at r u the relativistic-particle 
energy density dominates the magnetic energy density by 
roughly two orders of magnitude. The increasing particle 
dominance of Konigl jets with jet radius seems to be a com- 
mon feature of these models when n and m arc determined 
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the particle energy flux is of order 10 46 (1 + k) ergs s . 
This is of the same order as the particle kinetic energy 
fluxes found by Celotti & Fabian (1993), for a sample of 
sources that included 3C 279. The particle energy flux is 
about an order of magnitude higher than the bolometric 
radiative luminosity of 3C 279's jet (Hartman et al. 1996) 
after correction for beaming, which falls within the range 
of Lkin/Lrad found by Celotti & Fabian (1993). An en- 
ergy flux of 10 46 (1 + k) ergs s _1 is equivalent to an energy 
injection rate of order 0.1(1 + fc)M Q yr _1 , or a mass ac- 
cretion rate of order 0.1(1 + k)/r)Mo yr _1 , where 77 is the 
efficiency of conversion of mass to kinetic energy. 

We can compare the values of a, m, and n found here 
for 3C 279 with values of these parameters found for other 
sources where the Konigl model has been applied. The 
source apart from 3C 279 with the most constraints from 
spectral and VLBI data is 3C 345. Unwin et al. (1994) 
found a = —0.6, m = 1.5, and n — 1.4 for 3C 345 in mid- 
1990. Unwin et al. (1997) found a = -0.6, m = 1.9, and 
n = 1.7 for 3C 345 in mid-1992, and at this epoch they 
found that the Konigl-jet core was not the dominant X- 
ray emitter in the source. Hutter & Mufson (1986) found 
1.1 < m < 1.6 and 1.1 < n < 1.6 with an assumed a 
of —0.5 in their application of the Konigl model to three 
nearby BL Lac objects. 

4.2.3. Need for an additional homogeneous component 

The spectral index in the Konigl model fit for 3C 279 
that has the poorest observational constraints is the syn- 
chrotron index below the turnover frequency, a s \. This is 
because the extended jet emits a significant fraction of the 
flux at low radio frequencies, so that what we get from the 
single-dish spectra shown in Figure 6 is the spectrum of 
the core plus jet, when what we want is the spectrum of 
the core alone. The observed a s \ thus provides only a limit 
to the actual a s ±, with high- resolution VLBI at low fre- 
quencies being needed to accurately measure a s \. There 
are indications from VSOP observations at 1.6 and 5 GHz 
(Pincr et al. 2000) that, at least at that epoch, the spectral 
index of the VLBI core was much more inverted than the 
values of a s \ quoted in Table ^ If confirmed by further 
VSOP data on 3C 279 (Edwards et al. in preparation), 
this would rule out a pure Konigl jet model for the 3C 279 
core at these epochs. Because the Konigl model was cre- 
ated in part to explain the flat radio spectra of quasars, 
very inverted values of a s \ (less than about —1.0) create 
contradictions in the model, such as a synchrotron break 
that goes the 'wrong way' (a S 3 > a S 2). This situation 
could be rectified by adding a homogeneous component (a 
newly emerging 'blob') that is blended with the VLBI core 
on the VLBI images. 

It seems likely that, in the general case, the spectrum of 
the VLBI core is a blend of an inhomogeneous component 
like a Konigl jet and one or more shocks moving along 
the jet. Because these would all be merged on the VLBI 
images, it would be impossible to determine an inverse- 
Compton Doppler factor in this case, because the crucial 
observational constraint provided by the VLBI size would 
be missing. The degree to which the "quiescent" emission 
from 3C 279 can be represented by a single inhomogeneous 
component will determine the reliability of the Doppler 
factors computed from the conical-jet geometry. 



4.2.4. Speed and Orientation of the Jet 

The Lorentz factor 7 and angle to the line-of-sight 9 of 
the jet can be calculated if 6 and the apparent speed 8 app 
arc known: 

01 m + 5 2 + l , 
7 = app — (10) 



and 



= arctan 



26 

2/3, 



app 



B 2 +5 2 



1 



(11) 



In their similar work on 3C 345, Unwin et al. (1997) cal- 
culate the jet speed and angle to the line-of-sight at sev- 
eral radii along the jet. This was possible for 3C 345 be- 
cause the jet component C7 was the only good candidate 
for producing the X-ray emission, so the inverse-Compton 
Doppler factor measured for C7 (assuming C7 produced 
100% of the X-rays) could be combined with the apparent 
speed measured for C7 at different points along the jet to 
produce a plot of 7 and 9 vs. r (see Figure 4 of Unwin et 
al. 1997). 

The situation for 3C 279 is not so straightforward. As 
discussed in S 13.31 the core, inner-jet components, and C4 
all probably contribute a non-negligible fraction of the X- 
ray emission. A unique solution for 6 for each component 
can only be obtained by knowing a priori what this frac- 
tion is. Table gives Lorentz factors and angles to the 
line-of-sight for the conical-jet core under the assumption 
that the core produces all of the X-rays. Average val- 
ues of 7 and 9 obtained in this fashion are 7 = 6 and 
9 = 6°. Better estimates of 7 and 9 can be obtained by 
using instead the Doppler factor that reproduces the ob- 
served X-ray emission when all components are considered 
(see Table EJ). This provides a reasonable estimate for 7 
and 9 for the component listed as producing the majority 
of the X-ray emission at that epoch. These values of 7 and 
9 are listed in the final two columns of Table At epochs 
where the core or an inner-jet component was the domi- 
nant X-ray producer we used the average apparent speed 
of the inner jet of 5.2c. At the 4 epochs where C4 was the 
dominant X-ray producer we used the apparent speed of 
C4, or 7.5c (Paper I). In this fashion, we obtained average 
Lorentz factors and angles to the line-of-sight for the core 
and inner-jet region (r < 1 mas) of 7 = 8 and = 4°, and 
an average Lorentz factor and angle to the line-of-sight for 
C4 (at r w 3 mas) of 7 = 13 and 9 = 2°. These numbers 
apply only for a specific portion of C4's curved trajectory 
(see below), along this portion of its trajectory C4 is faster, 
in intrinsic speed and apparent speed, and its path is closer 
to the line-of-sight than the inner jet components. 

Whether this speed and angle are characteristic of C4's 
location in the jet or are unique to C4 is unknown, be- 
cause all components other than C4 faded by the time they 
reached 1 mas from the core during the period of our mon- 
itoring (see Paper I for a discussion). The quoted speed 
and angle for C4 apply for times when C4 was contribut- 
ing significantly to the X-ray emission, at the beginning 
and the end of the observed time range. The Doppler fac- 
tor of C4 is evidently time-variable and increased around 
1997, coincident with its brightening on the component 
light curves (Paper I). In Paper I we presented a detailed 
analysis of the kinematics of C4, based on its curved tra- 
jectory in the VLBI images. For that analysis we assumed 
7 = 13 (slightly higher than the minimum 7 required for 
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the maximum apparent speed along C4's curved path). 
We confirm this choice of 7 in this paper, and also confirm 
the other results from the kinematic analysis in Paper I: 
that the angle to the line-of-sight of C4 is about 2°, and 
that C4's Doppler factor was highest at the beginning and 
end of the observed time range (see Figure 8 of Paper I). 

5. CONCLUSIONS 

In this paper we have calculated, as accurately as is pos- 
sible, inverse-Compton Doppler factors for 3C 279. These 
Doppler factors were then used to compute the speed and 
orientation of the parsec-scale jet. Calculation of inverse- 
Compton Doppler factors is a notoriously inaccurate busi- 
ness. Nevertheless, given the large amount of multiwave- 
length spectral data recently published by HOI, and the 
large amount of VLBI data recently published by us (Pa- 
per I), 3C 279 seems to be the best object for constraining 
the Doppler factor by this method. Progress on this calcu- 
lation for other sources is hindered mainly by the reliance 
on the turnover frequency of the core and jet components, 
which is a critical parameter that is relatively poorly con- 
strained observationally, even for this well-observed source 
(see Figure 4) . Nearly simultaneous VLBI observations at 
as many frequencies as possible (and as near as possible to 
the same resolution) are required before attempting such 
a calculation. 

Major conclusions from this work are: 

1. The VLBI morphology and multiwavelength data 
cannot be adequately explained by either an inho- 
mogeneous jet (§ I3.2[l or a homogeneous sphere (or 
spheres) geometry (§ 13.1(1 alone. We have used a 
combined model (§ 13. 3[) . where we assume an inho- 
mogeneous conical-jet geometry for the VLBI core 
and a homogeneous sphere geometry for the VLBI 
components. In this combined model, we require the 
Doppler factor to be the same for the jet and spheres; 
the spheres can therefore be regarded as approxima- 
tions to dense clumps propagating along with the 



bulk jet. By applying this method, we obtain an 
average speed and angle to the line-of-sight for the 
core and inner-jet region (r < 1 mas) of v = 0.992c 
(7 = 8) and 6 = 4°, and an average speed and 
angle to the line-of-sight for C4 (at r « 3 mas) of 
v = 0.997c (7 = 13) and 6 = 2°. 

2. When the core is modeled as a conical jet, then on 
average the core produces about half of the X-rays, 
with the other half being split about evenly between 
C4 and the brightest inner-jet component. This re- 
sult differs from that for 3C 345 found by Unwin et 
al. (1997), who could not match the Konigl model 
to 3C 345 's X-ray emission, and concluded the core 
is not the dominant X-ray emitter in 3C 345. 

3. The jet is particle dominated at most radii that 
produce significant observed emission in the Konigl 
model. This result was also found by Unwin et 
al. (1994) for 3C 345. At the inner radius of the 
Konigl jet the magnetic field is of order 0.1 G and 
the relativistic-particle number density is of order 10 
cm~ 3 . The kinetic energy flux in the jet is of order 
10 46 (l + fc) ergs sec -1 , where k is the ratio of proton 
to electron energy, which implies a mass accretion 
rate of order 0.1(1 + k)/r)M & yr _1 , where 77 is the 
efficiency of conversion of mass to kinetic energy. 
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